Self-intersecting Regge trajectories in multi-channel scattering 
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We present a simple direct method for calculating Regge trajectories for a multichannel scattering 
problem. The approach is applied to the case of two coupled Thomas-Fermi type potentials, used 
as a crude model for electron-atom scattering below the second excitation threshold. It is shown 
that non-adiabatic interaction may cause formation of loops in Regge trajectories. The accuracy 
of the method is tested by evaluating resonance contributions to elastic and inelastic integral cross 
sections. 
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I. INTRODUCTION 

There has long been interest in resonance effects which 
arise when collision partners form a long-lived intermedi- 
ate complex pj . Recently, this interest was reinforced by 
experimental progress in cold atomic and molecular colli- 
sions (see, for example, [2]). An isolated resonance can be 
associated with a pole of the scattering matrix either in 
the complex energy (CE) or the complex angular momen- 
tum (CAM) plane. Two types of poles are closely related 
and contain essentially the same amount of information. 
However, since observables of interest such as integral 
(ICS) and differential (DCS) cross sections are given by 
partial wave sums (PWS) over total angular momentum 
J, the CAM (Regge) poles prove to be more convenient 
for identifying and quantifying resonance effects. Trans- 
forming a PWS into a sum of integrals, e.g., by means 
of the Poisson sum formula, and then evaluating contri- 
butions from Regge poles often allows to account for the 
resonance pattern observed in a DCS |3]-[T3j or an ICS 

m-m- 

Applications of the CAM approach range from elas- 
tic collisions of atoms with protons [M] and electrons 
[T5] , [12] to atom-diatom chemical reactions [13] , [TT] . Nu- 
merous techniques have been proposed for determination 
of Regge pole positions and residues in single channel 
potential scattering (for a review see [6]), among others 
direct solution of the Schroedinger equation [12] (see also 
[15)). or of the corresponding non-linear Milne equation 
[201 [5T] for complex values of J . For realistic reactive 
systems the number of open channel is large, root search 
in the complex J plane is not viable, and one has to re- 
sort to Pade' reconstruction of S"- matrix elements [H]. 
This leaves a class of systems with a relatively few chan- 
nels, for which many of the single channel techniques do 
not work, and yet one wishes to avoid the use of Pade' 
approximants. 

The purpose of the present paper is to propose a direct 



method for calculating CAM poles positions and residues 
for such systems. These include, among others, inelas- 
tic and reactive systems at low energies, Feshbach reso- 
nances, collisions involving two- level atoms [22], and spin 
flip scattering 23 . We will also look for the evidence of 
non-adiabatic effects in the behaviour of Regge trajec- 
tories. For recent efforts in this direction we refer the 
reader to Ref.[21], where the amplitude-phase method of 
Refs. |201 has been extended to Dirac electrons. The 
rest of the paper is organised as follows. In Sect. II 
we give a brief description of the method which gener- 
alises the approach of Ref. pTS] to a multichannel case. 
In Sect. Ill we consider a two-channel case designed to 
mimic electron-atom scattering below the second excita- 
tion threshold. Section IV considers the same problem in 
the adiabatic approximation. In Sect. V we use the ob- 
tained pole positions and residues to evaluate resonance 
contributions to elastic and inelastic integral cross sec- 
tions. Section VI contains our conclusions. 



II. DIRECT CALCULATION OF REGGE POLE 
POSITIONS AND RESIDUES 

Consider a time-independent scattering problem de- 
scribed by N coupled radial equations [we set to unity 
the particle's mass, /i = 1, and choose h = 1, thus con- 
verting to atomic units (a.u,)], 

{[-dl/2 + J{ J + l)/2r2 - E]i + V{r)}^{r) = (1) 

where E is the energy, J is the total angular momentum, 
/ is the unit matrix, and V{r) is an iV x hermitian 
potential matrix, such that 

Umr^ooV{r) ^ diag{Vi,V2, ...,Vn^uQ). (2) 

We will assume that the constant values Vn are arranged 
in such a way that Vi > V2 > ... > Vm-i, and re- 
quire that a solution of Eq.([T]), given by a complex vector 
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^(r) = [i'l, ^2, ^ni^ , is regular at the origin, 

;imr_,o*„(r) = 0, n = 1, 2, A^. (3) 

Furthermore, we are interested in scattering solutions 
^„j, TO = 1,2,...A'^, which for large r's contain an in- 
coming wave in only the m-th channel. Assuming that 
the potential V{r) reaches its asymptotic form ([2| suffi- 
ciently rapidly, as r oo for the channel wavefonctions 
we have 

« i7:Kr/2y^^[S„mH'^\knr) + (4) 
SnmiE,X)H['\knr)l n,TO=l,2...iV 

where h[^\z) and H^^\z) are the Hankel functions of 
the first and second kind, respectively, A = J -I- 1/2, and 
the asymptotic wave vector A:„ is given by 



1,2, ...TV. 



(5) 



For now we will assume that all channels are open, 
E > Vi so that all fc„ are real valued. For each real 
value of the energy E we wish to find complex value (s) 
A = X{E) , such that the S'-matrix elements Snm would 
diverge, Snm{E,X{E)) — oo. Thus, the asymptotic of 
the corresponding solution of Eq.Q (Regge state) will 
contain only outgoing waves generated by the emissive 
complex centrifugal potential [X{E)'^ — l/4]/2r^. 

In order to obtain Regge trajectory(s) X{E) we inte- 
grate Eq.([T]) for an arbitrary complex value of the an- 
gular momentum J (A) sufficiently far into the asymp- 
totic region, evaluate the S'-matrix elements, and repeat 
the procedure until a (Regge) pole of the S'-matrix is 
found [T5l HH] . With possible applications of the theory 
to electron-atom collisions in mind, we consider a poten- 
tial matrix which has a Coulomb singularity at the origin, 
so that its Taylor expansion takes the form 



V{r) 



3=0 



(6) 



As in the one-channel case [HKin] the singularity of V{r) 
prevents imposing the boundary condition ^ directly at 
r = 0. Following [TH] we represent the solution of Eq. 
(l)as a power series, 



(7) 



where Uj, j = 0,1,..., are constant vector coefficients 
satisfying recursion relations 

oo C30 

E ^j+uu + + 2{j + 1) E "i+iO' + (8) 



3=0 
oo j 



with Ki = Vi + EISi^i, whose explicit solution reads 
u,+i = -[2(J+ l)(j + 1) + j(j + 1)]"' > ^Kiu,.i, J > 0. 



3 

1=0 



The recursion scheme is initialised by specifying the so 
far undefined initial vector Uq. This can be chosen in 
N different ways, e.g., (the last subscript indicates the 
component of the vector u™) 

{ul^)n=Snm n, TO = 1 , 2, . . . , iV (10) 

to yield N linearly independent solutions, whose values 
at some vq for a suitably chosen jmax ■ 



*m(?'o) 



E"" 

3=0 



J+3 + 1 
'O 1 



(11) 



*™M= V(J + j + i)u 



3=0 



m J+j 
■] 'O ' 



TO = 1,2, ..N, 



provide N sets of initial conditions for Eq.([T]). This can 
now be integrated numerically (a NAG integrator 25 is 
used in this work) to a sufficiently large r where 



(*,„)„« (^fc„r/2)i/2[S 



TT 

mn X 



n, TO 



(2) 



{knr) 
1,2.. 



(12) 



and S~{E,X) and S'^{E,X) are constant matrices, to be 
determined numerically. The physical scattering states 
Q are linear combinations of 



(13) 



with Akm chosen so that the coefficients multiplying the 
Hankel functions of the second kind (incoming waves) 
add up to Smn i-G., A = (S^)^^. As a result, for the 
S-matrix we have 



^nn' — [{S ) 



(14) 



It is readily seen that the S-matrix elements diverge if 
and only if S" is singular, so that the condition for a 
Regge pole at a (real) energy E reads 



A(£:,A) = detS" 



0. 



(15) 



Starting with a reasonable initial guess for A and recal- 
culating the l.h.s. of Eq.(15) in each step, one can use 



a standard routine for finding zeroes of A{E, A) (a NAG 
root finder |26| is used in this work) to determine the 
accurate pole position A(£'). The residues 

p„„'(£;) = lim eS„„'(-B, A + e). (16) 

are readily obtained by integrating Eq.( T]) for a value of 
A close to X{E) and taking the limit (16 1. Finally, in the 



case some of the channels are closed, with corresponding 
fc„'s in Eq.([5| purely imaginary, equation (15) applies, 
and ensures that the Regge state does not have com- 
ponents which grow exponentially as r cxd. After this 
brief summary, in the next Sections we apply the method 
to a model two-channel {N = 2) problem. 



(9) 
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III. THE MODEL: TWO COUPLED 
THOMAS-FERMI TYPE POTENTIALS 

Next we consider a two-channel scattering problem 
with a potential matrix defined by 



1^22 (r) 



Z 



r(r + a)(r^ + b) 



and 



Vnir) ^ V22{r) + AV, 



Vi2{r) = = acxp[-(r - nf/Ar^ 



(17) 



(18) 



(19) 



The potential shown in Fig.l can be seen as a crude model 
for an inelastic electron-atom collision below the thresh- 
old of the second inelastic channel. The diagonal terms 
Vii and V22 are two similar Thomas-Fermi type poten- 
tials, representing the interaction between an electron 
and an atom in the first excited and the ground state, 
respectively. The constant AV is the excitation energy, 
and the interaction between the channels occurs in the 
outer layer of the atom, r,; — Ar ^ r < ri + Ar. Regge 




r(a.u.) 



FIG. 



1. (Color online) Elements of the potential matrix 
vs. r (solid). Also shown by dashed lines are the adiabatic 
potential curves 1 20 1 . 



trajectories (curves ImX vs. ReX) are shown in Fig. 2 for 
arbitrary Z = 54, a = 0.0125 a.u., h = 1.5874 a.u. and 
AV = 0.3 a.u. For uncoupled channels (V12 — V21 — 0), 
there are two Regge trajectories, Xi{E) and \ii{E). 
Since the two potentials only differ by a constant shift 
Ay, \i{E) = A//(£'-)- Ay), and the two Regge trajecto- 
ries in Fig. 2 coincide, as shown by the dot-dashed line. 
Interaction between the channels removes the degeneracy 
and yields two distinct trajectories as shown in Fig. 2 for 
a = 1.5 a.u., Ar = 1 a.u., and = 2.4 a.u. (sohd). The 
trajectory labelled (/) exhibits a type of behaviour often 
seen in potential scattering [141 EH HI] ■ As the energy 
increases, so does Xi{E), and the trajectory curves away 
from the real A-axis in a smooth manner. The second 
trajectory, labelled (//) leaves the real axis much more 




FIG. 2. (Color online) Regge trajectories ImJ vs. ReJ 
for the potential in Fig. 1 (solid) shown for E > AV. Also 
shown are Regge trajectories in the adiabatic approximation 
(201 (dashed) and the two coinciding Regge trajectories for 
''12 = V21 = (dot-dashed) 



rapidly, and shortly thereafter intersects itself, describ- 
ing a loop in the first quadrant of the CAM plane. This 
behaviour, to our knowledge not observed in one-channel 
scattering problems, is one of our central results. In the 
next Section we demonstrate it to be a consequence of 
non-adiabatic effects. 



IV. ADIABATIC CORRESPONDENCE 

Further insight can be gained by considering the two 
Regge trajectories in the adiabatic approximation. Di- 
agonalising the potential matrix V{r) for each value of r 
yields two adiabatic curves, 



Vi,ii{r) 



V1+V2 ^ ViVl - ^2)2 + 4^1 



± 



12 



(20) 



shown in Fig.l by dashed lines. The first (/) curve ac- 
quires an additional barrier, while the second one (//) 
has an additional well, both roughly proportional to 
±Vi2(r). Neglecting non-adiabatic coupling, i.e., replac- 
ing in Eq.Q 



V{r)~^diag[Vi{r),Vii{r)l 



(21) 



yields two uncoupled equations and two adiabatic Regge 
trajectories shown in Fig. 2 by dashed lines. The first 
adiabatic trajectory is close to the exact trajectory (/) 
in Fig. 2, whose imaginary part grows slowly with the 
energy, since the corresponding metastable state is sta- 
bilised by the effective barrier shown in Fig.l. The second 
exact trajectory (//) corresponds to the adiabatic trajec- 
tory for the barrierless potential V//(r). Comparing these 
two Regge trajectories suggests that self-intersection of 
the exact curve in Fig. 2 is caused by the non-adiabatic 
transitions not taken into account by the approximation 



(20l-(21) 
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FIG. 3. (Color online) Residue trajectories, Impij vs. Repij 
for (a) i=l, j=l; (b) i=l, j=2; (c) i=2, j=l; (d) i=2, j=2, for 
the trajectory (/) in Fig. 2 



V. INTEGRAL CROSS SECTIONS. 

One important application of the CAM theory is in 
the identification and quantitative analysis of resonance 
patterns which occur in elastic, inelastic and reactive in- 
tegral cross sections (ICS) [I4j-[T8]. A resonance is likely 
to affect an ICS at an energy E > for which the corre- 
sponding Regge trajectory approaches a real integer value 
of the angular momentum |14j.|15|. (The requirement is 
readily understood if one recalls that at negative ener- 
gies the condition for a true bound state is that J take 
a 'physical' integer value.) This condition is satisfied for 
the trajectory (I) in Fig. 2 which approaches J = 2 for 
E ^ 1.46 a.u. Next we use this example in order to check 
the accuracy of the pole positions and residues obtained 
in Sect. III. 

The four ICSs are given by the partial wave sums 



<Jnn'{E) = -5- V(2J-f 1)|(5„„- 



(22) 



j=o 



n,n' = 1,2 

where fc„/ is the wave vector in the incoming channel, and 



Sij is the Kronecker delta. The PWS (22 1 can be sepa- 



rated into the resonance and background contributions, 

a„n' (E) = (E) + Bnn' (E), (23) 

where the resonance term is given by the MuUhoUand 
formula [H], [TT], [IB] 



1 -I- exp(-27r?;A) 



(24) 



A = J -f- 1/2 and a * denotes complex conjugation. 

Exact composition of the background term -B„„' (E) is 
described elsewhere [TH [T7] , its essential property being 
smooth behaviour in the region where an ICS is affected 
by the resonance. Figure 3 shows the residues trajectories 
(curves Impij vs. Repij first introduced in Ref. [3]) eval- 
uated for the Regge trajectory (I) in Fig. 2. The full ICSs 
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FIG. 4. (Color online) State-to-state integral cross sections 
atj vs. E (thick solid) for (a) i=l, j = l; (b) i=l, j=2; (c) i=2, 
j = l; (d) i=2, j=2. Also shown are the resonance contributions 
alj" (solid) and the background terms Bij for the trajectory 
(I) in Fig. 2. 



given by the PWS ( 22 ) , the resonance term ( 24 ) , and the 
background term obtained as the difference between the 
two, are shown in Fig. 4. The resonance terms cr'^^'^E) 
account for most of the resonance structure at E k, 1.46 
a.u. We note that the trajectory (I) in Fig. 2 originates, 
at low angular momenta, from a bound rather than a 
metastable state of the two coupled wells. Thus, as in 
the case of proton impact on neutral atoms |14| and also 
electron-atoms collisions [16], one can expect Eqs.(23)- 



( 24 1 to provide an efficient separation of the resonance 
contribution and to probe important physics. Discussion 
of the distinction between two types of trajectories and 
a modification of the MulhoUand formula can be found 
in Ref. . Extension of the approach of Ref. [T5] to a 
multichannel case will be given elsewhere. 



VI. CONCLUSIONS AND DISCUSSION 

In conclusion, we advocate a direct method for cal- 
culating Regge pole positions and residues, suitable for 
systems with a relatively small number of channels. The 
method is applied to a simple model designed to mimic 
electron-atom scattering at energies between the first and 
the second excitation thresholds. It is shown that inter- 
channel coupling splits degenerate Regge trajectories into 
ones approximately corresponding to the two adiabatic 
potentials. Beyond the adiabatic approximation, non- 
adiabatic effects are seen to be responsible for self inter- 
section of the trajectory (//) shown in Fig. 2. The effect 
of loop formation has not, to our knowledge, been ob- 
served in single channel scattering. 

Finally, the simple model developed here can be im- 
proved e.g., by a careful choice of the potential matrix or 
by including, if necessary, additional J-dependent terms 
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in Eq.Q. The possibility to use the method for a more 
accurate description of inelastic electron-atom scattering 
will be discussed in our future work. Suffice it to say that 
this development promises a powerful approach to low- 
energy scattering with the possibility to probe Regge res- 
onances and the Feshbach resonances occurring in Bose- 
Einstein condesates . 
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